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ABSTRACT 

How and when the mass distribution of stars in the Galaxy is set is one of the main issues 
of modern astronomy. Here we present a statistical study of mass and density distributions of 
infrared dark clouds (IRDCs) and fragments within them. These regions are pristine molecular 
gas structures and progenitors of stars and so provide insights into the initial conditions of star 
formation. This study makes use of a IRDC catalogue (Peretto & Fuller 2009), the largest sam- 
ple of IRDC column density maps to date, containing a total of ^11,000 IRDCs with column 
densities exceeding Nh 2 = 1 X 10 22 cm~ 2 and over 50,000 single peaked IRDC fragments. The 
large number of objects constitutes an important strength of this study, allowing detailed analy- 
sis of the completeness of the sample and so statistically robust conclusions. Using a statistical 
approach to assigning distances to clouds, the mass and density distributions of the clouds and 
the fragments within them are constructed. The mass distributions show a steepening of the 
slope when switching from IRDCs to fragments, in agreement with previous results of similar 
structures. IRDCs and fragments are divided into unbound/bound objects by assuming Larson's 
relation and calculating their virial parameter. IRDCs are mostly gravitationally bound, while a 
significant fraction of the fragments are not. The density distribution of gravitationally unbound 
fragments shows a steep characteristic slope such as AN/ A\og(n) oc n~ 4 0±0 5 , rather indepen- 
dent of the range of fragment mass. However, the incompleteness limit at a number density of 
~ 10 3 cm -3 does not allow us to exclude a potential lognormal density distribution. In contrast, 
gravitationally bound fragments show a characteristic density peak at n ~ 10 4 cm~ 3 but the 
shape of the density distributions changes with the range of fragment masses. An explanation for 
this could be differential dynamical evolution of the fragment density with respect to their mass 
as more massive fragments contract more rapidly. The IRDC properties reported here provide 
a representative view of the density and mass structure of dense molecular clouds before and 
during the earliest stages of star formation. These should serve as constraints on any theoretical 
or numerical model to identify the physical processes involved in the formation and evolution of 
structure in molecular clouds. 

Subject headings: stars: formation — ISM: clouds 

1. Introduction bution of stars is determined is therefore essential 

in establishing a comprehensive picture of galac- 
While low mass stars dominate the mass of tic evolutiori; and star formation, throughout the 

galaxies, massive stars regulate their energy bud- Universe 

get. Understanding how and when the mass distri- „. „ . , , , , 

bmce stars form m molecular clouds the com- 
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Fig. 1. — Filamentary infrared dark clouds, (a) Spitzer 3 colour image of a region containing IRDCs (Blue: 
3.6/im ; Green: 8/im ; Red: 24/im ). The IRDCs are the long filamentary structures seen in extinction. 
On this figure the blue stars are foreground stars, the red/yellow stars are young stars currently forming in 
the IRDC. (6) H2 column density map constructed from the 8/im extinction seen in (a). While the outer 
contour delimits the boundary of each of the three IRDCs, the fragments are the substructures seen within 
each IRDCs. The 28 fragments identified in these clouds are marked by blue dashed ellipses. 



parison of the internal structure of the clouds 
and the initial mass function (IMF) of stars can 
provide insights on the processes responsible for 
the formation of stars. The mass distribution of 
molecular clouds, and cores within them have been 
extensively studied in the past twenty years. Until 
recently, it was believed that the mass distribution 
of CO clumps was described by AAco /A log M — 
M~ a with a = 0.7 ± 0.2 for the Milky Way 
(Kramer et al. 1998, Rosolowski 2004). The mass 
distribution of prestellar cores, the direct progen- 
itors of stars and stellar systems, as observed in 
dust continuum is much steeper, resembling the 
Salpeter IMF with a power law index of a — 1.35 
(Motte et al. 1998; Enoch et al. 2008). How- 
ever, recent papers questioned the impact of the 
source extraction scheme used to segment the data 
on the final mass distribution shape (Pineda et al. 



2009). Buckle et al. (2010) found a steeper mass 
distribution for small scale CO clumps. Also, in 
most cases, different tracers are required to trace 
different structures such as dense cores and molec- 
ular clumps, raising the question of detection bi- 
ases. Statistics is often a problem too, binning 
small number of objects introduce artifacts (Reid 
& Wilson 2006). Therefore some confusion exists 
on what is the real mass structure of molecular 
clouds. 

Another important physical aspect of molecular 
cloud structure is the probability density function 
(PDF) of the gas volume density. This quantity 
has received only little attention (e.g. Dring et al. 
1996 for HI; Smith & Scab 2009 for CO) but po- 
tentially contains crucial information on the pro- 
cesses at the origin of the density fluctuations. For 
instance, turbulence-driven fragmentation models 
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develop initial lognormal density fluctuations (e.g. 
Padoan et al. 1997), which could be the main 
driver of the lognormal part of the IMF (Chabrier 
2003). Studying the density distribution of frag- 
ments within molecular clouds could set important 
constraints on such models. 

To perform such studies, we decided to focus 
on a specific type of molecular clouds, i.e. infrared 
dark clouds (IRDCs). IRDCs are dense molecular 
clouds seen in silhouette against the bright emis- 
sion of the galactic plane (e.g. Perault et al. 1996; 
Teyssier et al. 2002; Rathborne et al. 2006; Simon 
et al. 2006). They are cold and only slightly pro- 
cessed by star formation activity, still containing 
the initial conditions of star formation. Peretto & 
Fuller (2009; hereafter PF09) recently constructed 
the column density maps of more than 11,000 of 
such IRDCs, the largest database of such struc- 
tures to date. This catalogue provides the oppor- 
tunity to probe molecular clouds in the Galaxy 
over a wide range of size scales and column den- 
sity at high angular resolution using the 8//m dust 
absorption and a new source extraction scheme. 

In section 2 of the present paper we discuss the 
dataset we used. In section 3 we describe and 
re-analyze previous results on the distance distri- 
bution of IRDCs, while in Section 4 we estimate 
completeness limits. Section 5 displays our main 
results on the size, density and mass distributions 
of IRDCs and fragments. Discussion is in Section 
6. And finally we summarize the main findings of 
this paper in Section 7. 

2. Data set 

The analysed IRDCs come from a new cata- 
logue of clouds identified in the Spitzer GLIMPSE 
data (PF09). IRDCs were defined as connected 
structures with column density peaks above 
JVh 2 = 2 x 10 22 cm~ 2 and boundaries defined 
by the contour at Nn 2 = 1 x 10 22 cm~ 2 . Sin- 
gle peaked structures lying within the IRDCs 
were identified as fragments. The boundary of 
a fragment being defined by the contour of the 
local minimum between a fragment and its closest 
neighbour, the same criterion used to define the 
leaves of the dendogram analysis of Rosolowsky et 
al. (2008). As column density peaks, these frag- 
ments are particularly important in the context of 
star formation since they are the likely birth place 
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Fig. 2. — Histogram of the IRDC peak column 
density. We can see a steady decrease down to 
Nn 2 ~ 1 x 10 23 cm -2 . Above this limit saturation 
does not allow us to probe the true peak column 
density. 

of the future generation of stars. The catalogue in- 
cludes opacity maps at 4" resolution and physical 
properties for over 11,000 IRDCs. Extracting the 
densest structures, a total of ~50,000 fragments 
have been catalogued within the full sample of 
clouds (PF09). Figure la shows a Spitzer three 
colour image of a region containing three filamen- 
tary IRDCs from the PF09 catalogue. Figure lb 
shows the column density map of these IRDCs 
and identifies the fragments within the clouds. 

2.1. IRDC saturation 

As discussed in PF09, some of the absorption 
towards IRDCs is saturated, meaning the infrared 
background is not strong enough in order to fully 
probe the internal structure of an IRDC. Based 
on photometric noise limitation and background 
strength PF09 estimated the fraction of saturated 
IRDCs to be 3%, corresponding to roughly 340 
IRDCs over the entire sample. 

An alternative estimate of the number of sat- 
urated clouds can be made from an inspection of 
the distribution of peak column density of IRDCs 
shown in Fig. [2] There appears to be a break in 
the distribution of peak 7Vh 2 at ~ 1 x 10 23 cm~ 2 , 
which likely reflects the effect of saturation. The 
fraction of IRDCs lying above this limit is 4%, very 
similar to the value estimated in PF09. However, 
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even for these saturated clouds only a small frac- 
tion of their area is above the saturation limit, only 
marginally affecting the averaged IRDC column 
density (and therefore any estimate of the cloud 
mass). But, the saturation has a much stronger 
effect on some fragments. For this reason, in the 
analysis presented here IRDCs containing satu- 
rated pixels are considered, but fragments with 
saturated pixels are excluded. 

Our estimated saturation limit is roughly twice 
as large as the one found by Vasyunina et al. 
(2009) from millimetre emission in their study of 
particularly high column density IRDCs. The dis- 
crepancy between the low absorption column den- 
sities Vasyunina et al. determined by assuming 
the minimum possible foreground emission, that 
due to the zodical light, and the high values they 
determined from millimeter dust continuum led 
Vasyunina et al. to derive a relatively low satura- 
tion limit. However, the majority of their clouds 
do not in fact appear saturated as considerable 
substructure can be seen in the 8/mi extinction 
maps. 

2.2. Column density and angular size dis- 
tributions 

This study aims to statistically analyze the den- 
sity and mass distributions of IRDCs and their 
fragments. To derive such quantities we first 
need to know the angular size and column density 
distributions as measured on the column density 
maps constructed by PF09. 

Figure [3] shows the distribution of the angular 
size and column density for the ~ 11,000 IRDCs 
and the ~ 50, 000 fragments identified within 
them. Fragments with saturated pixels have been 
excluded (see Sec. 2.1). In addition IRDCs which 
are not fragmented (~ 40% of the IRDC sample) 
have also been removed to maintain a clear def- 
inition of a fragment as a substructure within a 
cloud. However, in practice keeping these single 
peak clouds has little effect on the results. 

It is important to note that the column den- 
sities we plot here are the background substracted 
column densities, equivalent to the one obtained 
in the clipping option of the dendogram analy- 
sis of Rosolowsky et al. (2008). In the context 
of centrally concentrated structures, these column 
densities are the relevant ones when interested in 
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Fig. 3. — Distributions of angular radius R eq and 
average column densities over each structure for 
IRDC (top) and fragments (bottom). The esti- 
mates for the radius completeness limit are dis- 
cussed in Section HI 

the physical properties of the gas enclosed in a 
given radius. Figure |3] clearly shows that the dis- 
tributions are dominated by small structures of 
low column density. We can also clearly see the 
effect of incompleteness on the distributions with 
the decrease in the number of sources at low ra- 
dius/column density, responsible for the formation 
of artificial peaks. The incompleteness in the sam- 
ple and these distributions are discussed in Sec- 
tion n 

3. IRDCs distance distribution 

To calculate the density and mass of the clouds 
the distance of each IRDC is required, however 
this is not yet known for most of the 11,000 IRDCs. 
For this analysis we have therefore adopted a sta- 
tistical approach based on previous measurements 
of the distances to samples of IRDCs. 

Several studies have measured the distance dis- 
tribution of subsamples of IRDCs in both the 
1st and 4th quadrant of the Galactic plane (Si- 
mon et al. 2006; Jackson et al. 2008; Marshall 
et al. 2009). Both kinematic and dust extinc- 



4 




Jackson et al. (2008) 
4th quadrant 

Simon et al. (2006) 
1 st quadrant 



Marshall et al. (2009)_ 
4th quadrant 

1st quadrant 




B 8 



o 
o 



2 4 6 8 10 
IRDC heliocentric distance (kpc) 



-m- 



1 



□ [] 



[] □ 



-50 

Galactic longitude (degree) 



50 



Fig. 4. — (left): Distributions of distances of IRDCs (Simon et al. 2006, Jackson et al. 2008, Marshall 
et al. 2009). (right): Plot of the galactic longitude dependance of average IRDC distance determined 
from kinematics (top panel) and extinction (bottom panel). The horizontal dashed lines show the average 
distances over the indicated longitude range (Table Qj) . 



tion techniques have been used to infer these dis- 
tances, and although they lead to similar results, 
there are some differences (Fig. @|. The proper- 
ties of these distance distributions are summarised 
in Table [1] Using dust extinction, Marshall et 
al. (2009) found a centrally peaked, Gaussian- 
like distance distribution very similar for both the 
1st and 4th quadrant of the Galaxy (Fig. @] - bot- 
tom panel) . In a similar way, kinematic distance^] 
shown in Fig. 0] (top panel) show a good agreement 
between the 1st and 4th quadrant . Although in 
the 1st quadrant a tail at 5 kpc clearly emerges. 
However most significant difference between the 
extinction and kinematic distances is the position 
of the peak, being located at ~ 3 kpc in one case 
and at ~ 5 kpc in the other. Both techniques have 
their own biases and advantages, it is therefore 
difficult to favor one distance distribution over an- 
other. However a Gaussian distribution is a rather 
good approximation to the distance distribution in 
both quadrants. 

Fig. E] (right) shows the average distance to the 
IRDCs as a function of galactic longitude for the 



These kinematic distances have been recalculated by taking 
the CS(2-1) and 13 CO(1-0) velocities published in Jackson 
et al. (2008) and Simon et al. (2006), respectively, and 
using the Reid et al. (2009) revised galactic rotation model. 



sources with distances measured by these two tech- 
niques. Any systematic trend with longitude could 
introduce a bias in analysis adopting a statistical 
distribution for the distance of IRDCs. It is clear 
that there is very little variation of the IRDC dis- 
tances with respect to the galactic longitude. The 
only region where there may be such an effect is 
towards galactic centre seems, an area which is not 
covered in our IRDC Spitzer catalogue which only 
extends into I = ±10°. It is worth noting that dis- 
tance variations around the mean as a function of 
longitude are very similar for both methods, em- 
phasising that is it predominantly only the average 
distance which differs between the two methods. 

Another possible bias is with respect to the 
size of the IRDCs. The studies shown in Fig. [4] 
do not contain IRDCs as small as those in the 
Spitzer based sample and so it is possible that 
the small and large clouds have different distance 
distributions. However, recent observations with 
ATNF Mopra in CS J=l-0 (Peretto et al. in 
prep) of a square degree of the galactic plane 
(29.8° < I < 31.8°, -0.27° < b < 0.27°) can be 
used to investigate this possibility. This area cov- 
ers 196 Spitzer IRDCs in total and more than 80% 
of the clouds can be associated with CS emission. 

The distances of these clouds can be calculated 
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Sample 


Mean distance 
(kpc) 


FWHM 

(kpc) 


Dispersion, a 
(kpc) 


Marshall, 1st & 4th quad. 


4.7 


1.5 


1.2 


Jackson, 4th quad. 


3.0 


0.8 


1.2 


Simon, 1st quad. 


3.6 


2.5 


1.3 



Table 1: Properties of the distance distributions of 
Gaussian the FWHM would equal 2.35o\ 



t 1 1 1 1 1 1 1 1 1 r 
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Heliocentric distance (kpc) 

Fig. 5. — Heliocentric distance distribution of 
IRDCs observed in CS J=l-0 (Peretto et al. in 
prep). The 3 histograms correspond to 3 different 
range of IRDC size. There is no evidence of any 
statistical difference in IRDC distance as a func- 
tion of their size. 

using the Reid et al. (2009) galactic rotation 
model. Figure [5] shows the distance distribution 
of the Spitzer IRDCs detected in CS. In this fig- 
ure, the IRDCs have been divided int to three size 
ranges. The 80 smallest IRDCs, those with R eq 
(PF09) less than 15" have a mean distance and 
standard deviation of 5.2 kpc and 0.8 kpc. This 
is indistinguishable from the values of 5.3 kpc and 
0.7 kpc and 5.3 kpc and 0.8 kpc for the IRDCs 
in the next two size ranges, 15" < R eq < 30", 
R eq > 30", which contain 52 and 39 objects re- 
spectively. These distributions therefore show no 
indication that large and small IRDCs have differ- 
ent distributions of distance. 



's shown in Fig. If the distributions were truly 



4. Completeness limits 

4.1. Column Density and Mass 

The mass completeness limit for the IRDCs and 
fragments can be written as 

M c = Tr(R c d c ) 2 x <N H2 > c (1) 

where R c is the smallest radius above which the 
sample is complete, d c is the distance within which 
the majority of the sources occur, and < Nn 2 > c 
is the typical average column density of the struc- 
tures with a radius R c . Figure |4]deft) shows that 
about 95% of the IRDCs in that plot have dis- 
tances below 6 kpc and so we conservatively adopt 
d c = 6 kpc. 

Estimating R c is less straighforward. The com- 
pleteness limits of our survey are related to two 
parameters of the source extraction: iV™ p , the 
minimum column density amplitude of a source 
(which is related to sensitivity) from the boundary 
of a cloud to its peak; and the angular resolution, 
4" both for IRDCs and fragments. In order to in- 
vestigate how these contribute to the completeness 
limits we look at the distribution of average col- 
umn density of IRDCs for objects of a given range 
of sizes as plotted in Fig. [51 We then plot the col- 
umn density at the peak of these distributions as 
a function of cloud size. This is also done for the 
fragments. Figure [7] shows these plots. 

The plots show a similar structure for both the 
IRDCs and fragments. Up to some size, 55" for the 
IRDCs and 9" for the fragments, the peak of the 
column density distributions is constant. Above 
these values it increases with increasing size. This 
constant column density for small sizescales sug- 
gests that the sample is not fully probing the pop- 
ulations of objects at these sizescales. There are 
objects in these size ranges which have lower col- 
umn densities and are not sampled by the objects 
in the catalogue. These plots therefore show the 
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Fig. 6. — Peak column density distributions for different IRDC (left) and fagment (right) sizes. For both 
structures we show examples for which we are sensitivity limited and for which we are not. The location of 
the peaks of these distributions as a function of size of the clouds and fragments are used to construct Fig. [7] 
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Fig. 7. — Maximum peak column density for all IRDCs (left) and fragments (right) falling in a given bin 
of radius. We can see that this maximum differs from the sensitivity limit N-^ when reaching a certain 
radius. This radius is taken as being the completeness radius. 



size limit completeness of the catalogue, for both 
IRDCs, R c = 55", and fragments, R c = 9". 

Adopting these sizes, the average column den- 
sity of clouds/fragments below these sizes gives 
average column density completeness limits of < 
N H2 >f' ag = 2.5 x 10 21 cm" 2 and < N H2 >J. RDC = 
5 x 10 21 cm -2 . Using Eq. [TJ these values give the 
mass completeness of the catalogue as M* = 
800 M Q and M c fra « = 9 M Q . 

4.2. Density 

The density distribution of IRDCs is difficult 
to interpret since the clouds are defined based on 
a column density threshold. Therefore we confine 
our discussion of the density distributions to the 
fragments. Both sensitivity and angular resolu- 



tion are important limiting factors in the context 
of density completeness of the sample: both com- 
pact low-mass fragments and large, diffuse, high- 
mass fragments could remain undetected. For any 
fragment mass, M\i m there is a minimum radius 
for which the peak column density of the fragment 
becomes higher than the threshold for identifying 
fragments (3 x 10 21 cm -2 ). If this minimum radius 
is larger than the angular resolution then such a 
fragment is detected. Therefore, we can define a 
density completeness limit for all fragments more 
massive than Mu m . The minimum radius, i? m i n , 
and the corresponding maximum density, p l ° w , are 
given by 

p _ / M lim / 9 x 

A mm — -l / j\r am P ^ V > 
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where fi is the mean mass per molecule and we 
adopt a column density < N^ p > which is the 
average column density of fragments of mass My lm 
(which by definition have peak column densities 
greater than the threshold to be identified as frag- 
ments). At a distance of 6 kpc, for fragment mass 
of 2/8/32 M and < N^ p >= 1.5 X 10 21 cur 2 
we get R min = 5/10/20" and p[ ow = 2.4/1.2/0.6 x 
10 3 cm -3 , respectively. 

For the same given mass there is also an up- 
per density limit which corresponds to the point 
where the size of the fragment becomes smaller 
than the resolution of the observations. This 
provides an upper limit on the density, = 
0.75 x Muraj '(■nR^. es ). For the 3 masses discussed 
before we get p%*> = 0.4/1.7/6.8 x 10 5 cur 3 . 

5. Size, mass and density structure of 
IRDCs 

5.1. Physical size distribution 

To calculate the density and mass of the IRDCs 
and fragments requires a distance for each ob- 
ject. Two different approaches to statistically at- 
tribute a particular distance to a particular cloud 
have been adopted. The first is to simply as- 
sign a unique distance to all clouds. Doing this, 
the physical size distribution of IRDCs and frag- 
ments will be exactly the same as the angular size 
distribution. Given the well peaked distribution 
of distances for clouds with measured distances 
(Fig. U]), this should be a reasonable first approx- 
imation. However, a more sophisticated approach 
is to make use of the distribution of distances 
(rather than just its peak position). To do this 
we adopt a distance distribution for the IRDCs 
and then randomly assign a distance drawn from 
this distribution to each cloud. Doing this for the 
whole sample of clouds repeatedly provides a sta- 
tistical sampling of the distance distribution. The 
final physical size distribution is the convolution 
of the true physical size distribution by the chosen 
distance distribution. However this does not have 
a crucial impact on the interpretation of the phys- 
ical size distribution if the dispersion of the dis- 
tance distribution is much smaller than the angu- 
lar size distribution. This is clearly the case since 



the angular sizes, both for IRDCs and fragments, 
extends over 2 order of magnitudes while IRDC 
distances span only over a factor of 3 at most. In 
other words, the dispersion in distance has rela- 
tively little effect on the final physical size distri- 
bution (see Appendix A). To assign distances to 
the clouds using this sampling technique we adopt 
a Gaussian distribution of distances with a peak 
at 4 kpc and a dispersion of 1 kpc, consistent with 
observed distance distributions (Fig. 0]). 

5.2. Mass distributions 

Mass distributions of molecular cloud struc- 
tures have been extensively studied in the past, 
therefore they represent a good point of compari- 
son for this current study. We defined mass as: 

M = 7rR 2 cq < N H2 > (4) 

where < Na 2 > is the average column density 
across the IRDC or fragment and R eq its equiv- 
alent radius (PF09). 

Figure [5] shows the mass distributions for 
IRDCs and fragments calculated adopting a sin- 
gle distance of 4kpc (filled square symbols) and 
for randomly attributed distances as described in 
Section [Q1 The shaded band on the figures shows 
the range (3 times the dispersion) spanned by the 
100 different distance realizations and the open 
triangles the mean for the different realizations. 
The completeness limits are shown by the dashed 
lines. For comparison the power-law slopes of the 
CO clump mass function (slope = —0.7) and the 
Salpeter mass function (slope = —1.35) are also 
shown. Using the MPFITS IDL package (Mark- 
wardt 2009) we have fitted the two distributions 
above their respective completeness limits. For the 
IRDCs we find a linear function (in a log-log plot) 
provides a good fit with dNiBDc/d log M — M~ a 
with a = 0.85 ± 0.07. The mass distribution of 
fragments is better fitted by a lognormal function 
defined as 

^|2£ = Aexp(-pog(M) - log(M pcak )] 2 /2 ( x 2 ) 

(5) 

where A is a normalization constant, Af pca k is the 
peak mass of the distribution, and a is the dis- 
persion. However, since we do not map the peak, 
the precise parameters of the lognormal function 
fit are not well constrained, several provide ade- 
quate fits to the data points. The function showed 
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Fig. 8. — Mass distribution of infrared dark clouds (left) and fragments (right) for two different distance 
distributions: The filled square symbols correspond to the adoption of an unique distance of 4 kpc for each 
single cloud, while the open triangles and associated shaded area corresponds to adopting and sampling a 
Gaussian distance distribution (see Section 5.1). The vertical dashed lines show the incompleteness limits. 
The best-fit is linear for the IRDCs (red solid line) with AAi RD c/A log(M) = M~ a with a = 0.85 ± 0.07, 
while the best-fit for fragments (blue solid line) is a lognormal function. For comparison the slope of the 
mass distribution of CO molecular clouds and clumps and the Salpeter part of the IMF are also shown. 



in Fig. E] has A=4610, M peafe = 1.55 M , and 
a = 0.78. As argued above, the figure confirms 
that as a consequence of the nature of the dis- 
tance distribution, there is relatively little differ- 
ence in the derived mass distributions whether a 
single distance is adopted for the clouds or statis- 
tical approach is adopted. Also, the results of this 
analysis are also not strongly dependent on exact 
parameters of the assumed distance distribution 
as demonstrated in Appendix A. 

A number of previous studies have attempted to 
construct, with samples at least an order of mag- 
nitude smaller, the mass distributions of IRDCs 
(Simon et al. 2006; Marshall et al. 2009) and 
fragments within them (Rathborne et al. 2006; 
Ragan et al. 2009). Except for the Ragan et 
al. study, the mass distributions in these studies 
agree: the IRDC mass distribution is similar to 
that of CO clumps, while the distribution for the 
sub-structures are steeper, more like the Salpeter 
IMF. 

In their analysis of 11 IRDCs, Ragan et al. 



(2009) found that the mass distributions of what 
they called clumps, which correspond to fragments 
here, is quite flat, similar to the CO clump mass 
distribution, in contrast with the present study. 
However it is difficult to understand the Ragan et 
al. result as the radii and masses they quote for 
their clumps imply 8/im opacities over 10 times 
larger than the 8/xm opacities they quote. 

5.3. Density distribution 

The density distribution of fragments may pro- 
vide important insights on the physical process 
generate these structures. We define the number 
density of a fragment as 

n =< N H . 2 > /dt (6) 

where d t is the line of sight size of the fragments 
which is assumed to be twice the projected radius. 
Figure [5] shows the fragment density distributions 
for the following mass ranges: 2 < M < 8M Q , 
8 < M < 32M ra and M > 32M . For each 
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Fig. 9. — Distribution of the number density of fragments normalised to 100 cm~ 3 . (left): Fragments with 
a mass 2 < M < 8 M©. (middle): Fragments with mass 8 < M < 32 M©; (right): Fragments with a mass 
> 32 Mq. The dashed lines mark the density completeness limits, lower and upper for the left hand side 
panel and only lower for the 2 others (since the upper limits are off the plots). The square and triangle 
symbols and shaded area have the same meaning as for Fig. [8] 



range, the density completeness limit, the dashed 
line, is calculated for Mn m = 2/8/32 M©, respec- 
tively. The density distribution over the entire 
mass range is very similar to that for the lower 
mass range (left panel). The figures show that go- 
ing from low mass to high mass fragments, the 
distributions become flatter. Compared to low 
density fragments, there is a higher probability of 
finding high density fragments for high mass frag- 
ments. One of the main issue in interpreting such 
a plot in terms of the formation of the fragments 
is that the density of gravitationally bound frag- 
ments is increasing over the time as they evolve 
(i.e. contract), and therefore might pollute the 
initial density PDF of the parental IRDCs. In the 
next section we will discuss the impact of such ef- 
fects on the density distributions. 

6. Discussion: Turbulent vs gravitational 
dominated structures 

The mass distributions of IRDCs and fragments 
plotted in Fig. [5] clearly shows a steepening, from 
large structures to smaller fragments. While the 
mass distribution of IRDCs is similar to that of CO 
clumps, the fragment mass distribution has a slope 
at high masses which is reminiscent of the slope of 
the Salpeter IMF, although it is best fitted with 



a lognormal function. However, two biases could 
affect the shape of the high-mass end of the distri- 
bution and the interpretation that its slope is re- 
lated to the Salpeter IMF. The detailed structure 
of this part of the distribution may be particularly 
sensitive to the adopted Gaussian distance distri- 
bution. Also high-mass fragments might evolve 
more rapidly than their low-mas analogues and 
therefore be under represented in extinction ob- 
servations at 8/xm (c.f. Hatchell & Fuller 2008). 

To the first order, the mass distributions are in 
agreement with the theoretical work of Hennebelle 
& Chabrier (2008) who interpreted the transition 
from a flat mass distribution to a steeper one 
as the transition from turbulence-generated struc- 
tures to gravity dominated structures. In this 
context it is interesting to measure the gravita- 
tional binding of IRDCs and fragments. To inves- 
tigate this, we have used Larson's relation (Larson 
1981) to compute the kinetic support, and cal- 
culate the virial mass. Following Hennebelle & 
Chabrier (2008) we assume the effective velocity 
dispersion is given by 

ccff = {cl+QMV^d 2 ^) 1 ' 2 (7) 

where c s is the sound speed Vo is the normaliza- 
tion velocity of Larson's relation, rj is the power 
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Fig. 10. — Fraction of bound IRDCs and frag- 
ments as a function of their mass. The typical er- 
ror due to distance uncertainty is 0.1. The dashed 
lines show the incompleteness limits of both frag- 
ments and IRDCs. 

law index of Larson's relation, and d is the size 
of the structure. We can then use this to com- 
pute the corresponding virial mass M v ; r for every 
IRDC/fragment. Figure [TOl shows the fraction of 
IRDCs and fragments having a ratio Ma 2 /M vir > 
0.5, assuming c s = 0.2 km/s (T=10 K), V = 
1 km/s, and rj — 0.4. Above the completeness 
limit all the IRDCs appear gravitationally bound, 
as do the majority of the fragments. Of course 
large uncertainties exist on the use of Larson's re- 
lation and its normalization. However, different 
normalizations still give similar conclusions about 
the fraction of IRDCs and fragments which are 
bound. 

A consequence of the IRDCs being gravita- 
tionally bound is that the observed change of 
slope of the mass distributions shown in Fig. |8] 
does not represent the transition from turbulence 
dominated to gravity dominated structures: most 
of the IRDCs down to the completion limit are 
bound. However, as shown in Fig. [10] a signif- 
icant fraction of the fragments lying above the 
completeness limit are unbound. In other words, 
even if globally gravitationally bound, IRDCs may 
contain turbulence-generated over-densities which 
will probably disperse and not form stars. The 
physical properties of these unbound fragments 
are likely to be representative of the initial con- 
ditions of star formation within IRDCs. 

Using the previously defined ratio to separate 



bound fragments to unbound ones we constructed 
the density distributions for both type of frag- 
ments as shown in Fig. [TTJ The mass ranges are 
the same as in Fig. The unbound fragments 
have all very similar distributions, independent of 
their mass range. In particular, the high mass end 
of the distribution is well fit by the following re- 
lation AA/Alog(n) oc n~ 4 ' 0±0 ' 5 , the error bar 
arising from the uncertainties on Larson's 
relation. The location of the peak seems to move 
with the completeness limit and is therefore ques- 
tionable. From these plots we cannot exclude a 
possible lognormal distribution for unbound frag- 
ments, but the peak of such distribution would 
have to be below n~ 10 3 cm~ 3 . 

In contrast, the density distribution of the grav- 
itationally bound fragments show a well defined 
peak between n=10 3 cm -3 and 10 4 cm' 3 in each 
mass range and a shape which broadens to lower 
densities as the mass range increase. This could 
result from the higher mass fragments of a given 
density evolving to being bound more rapidly that 
lower mass fragments. 

7. Summary 

We used the largest sample of IRDC column 
density maps to date in order to better charac- 
terize the size, mass, and density structure of 
dense molecular clouds. The large number of 
objects, 11,000 IRDCs and 50,000 fragments, al- 
lows a detailed analysis of the completeness of 
the sample. Using a statistically attribute dis- 
tances to each IRDC, we have demonstrated that 
above the completeness limit the mass distribu- 
tion of the IRDCs are consistent with a power-law 
AJVi R Dc/Alog(M) = M -0,8 , where N mDC is the 
number of clouds. For the fragments the high mass 
end of the mass distribution shows a steeper slope, 
consistent with the slope of the Salpeter IMF, with 
the overall distribution well fitted by a lognormal 
function. 

Using Larson's law to estimate the linewidth 
of each IRDC and fragment, we have shown that 
above our completeness limit all the IRDCs and 
the majority of fragments are likely to be bound. 
This implies that the transition in the shape of 
the mass distribution does not reflect a transition 
from unbound to graviationally bound structures. 
Looking at the distribution of fragment density 
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Fig. 11. — Same as Fig. [HI the shaded areas are the same as in Fig- EI The red solid line represent the density 
distributions of gravitationally unbound fragments, while the blue dashed line represent the gravitationally 
bound fragments. 



shows that bound fragments dominate the high Patrick Hennebelle for some useful discussions. 

density (n > 10 4 cm -3 ) end of the distribution for 
all mass ranges, and dominate the whole distri- 
bution for the highest range of fragment masses. 
There is also a distinct broadening of the distri- 
bution with increasing fragment mass. This could 
be a result of the higher mass fragments evolv- 
ing to being bound more rapidly that lower mass 
fragments. The number of unbound fragments 
as a function of number density has the form 
AiVf/Alog(n) oc n~ 4 - 0±0 - 5 (where Nf is the num- 
ber of fragments) down to a density of ~ 10 3 cm~ 3 
where the completeness limit is reached. 

The absence of bright infrared sources embed- 
ded in IRDCs indicates that the mass distributions 
and density distributions as a function of mass and 
degree of gravitational binding derived here are 
representative of the initial conditions of star for- 
mation within dense molecular clouds. These re- 
sults should serve as constraints on theoretical and 
numerical models in order to identify and charac- 
terize the physical processes responsible for the 
formation and early fragmentation of molecular 
clouds. 
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A. Different distance distributions 



In order to investigate further the impact of choosing a given distance distribution on the calculated mass 
distributions we performed a series of tests with different assumptions. Using the same approach as described 
in Section 4 and 5, we used a Gaussian distribution for distances with a mean distance of 4 and 5 kpc with 
dispersions of 0.5, 1 and 2 kpc. We also performed a test assuming different distance distributions for the 
1st quadrant and 4th quadrant IRDCs. We show the resulting mass distributions in Fig. rT2"land[T3l We see 
that the shape of the mass distributions above the completeness limits for both the IRDCs and fragments, 
remain basically the same, with best fit parameters changing only of a few 0.1 dex. 
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Fig. 12. — Mass distribution of IRDCs. The different symbols have the same meaning as in Fig. [S] The 
four plots corresponds to four different assumptions on the IRDC distance distributions: (top-left) Gaussian 
distribution with a peak at 4kpc and a dispersion of 0.5kpc; (top right) Gaussian distribution with a peak 
at 4kpc and a dispersion of lkpc; (bottom-left) Gaussian distribution with a peak at 4kpc and a dispersion 
of 2kpc; (bottom-right) Gaussian distribution with a peak at 5kpc and a dispersion of 2kpc for the 1st 
quadrant clouds and a Gaussian distribution with a peak at 4kpc and a dispersion of lkpc for the 4th 
quadrant clouds. Note that we also performed tests with a peak distance at 5kpc for all IRDCs but these 
resulted in distributions very similar to the 4kpc results shown here, but shifting the distributions slightly 
to more massive objects. 
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Fig. 13. — Same as Fig. H"2l but for the fragments. 
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